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(54) DS-CDMA RECEIVER 

(57) A path search stage comprises smaller number 
of sliding correlators than the code length of a predeter- 
mined spreading code, calculates a correlation at a 
phase having a center of correlation peaks detected 
with a matched filter as a center timing, and causes 

RS , 1 ■ 



change in the calculation frequency according to the 
change in a delay profile. Thereby, the reduction in the 
circuit size and the electric power consumption may be 
realized. 
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M , .I- 

Description 
Technical Field 

5 [0001] The present invention relates to a receiver for a direct sequence code division multiple access (DS-CDMA) 
communication to transmit a signal obtained through multiplication of an information signal for transmission by a prede- 
termined spreading code at the transmitting part, and calculate a correlation between a received signal and the above 
spreading code and demodulate the above information signal at the receiving part. 

w Background Art 

[0002] Conventionally, a receiver of a DS-CDMA system, more particularly a receiver of a mobile station has been 
provided with a cell search stage to perform initial cell search, and a path search stage to detect a multipath timing. The 
cell search stage calculates a correlation between a received signal and a spreading code, and obtains a correlation 

?5 power to select a base station as a transmitting source. On the other hand, the path search stage calculates a delay 
profile by correlation operation between a received signal and a spreading code to detect a multipath timing (phase) 
based on the delay profile. FIG. 15 shows a format of a radio transmission frame of a received signal processed at the 
path search stage. The radio transmission frame of the signal comprises a plurality of slots, and each slot SLOTk has 
a pilot block PILOTk and a data block DATAk. The pilot block PILOTk includes a known signal, based on which fading 

20 compensation of the data block is performed. Thereafter, the delay profile of the data block is calculated. 

[0003] And data blocks and pilot blocks in FIG. 15. in some cases, maybe transmitted as two parallel channels, that 
is. maybe divided to two radio frames, one comprising only data blocks, and the other only pilot blocks. In the above 
cases, each user of the same cell performs fading compensation of the data blocks on the known signal for the channel 
comprising only pilot blocks, and the channel is called as a common pilot channel. 

25 [0004] The mobile station has been strongly required to have the minimum eleciric power consumption in order to 
secure the portability and long operation time. Then, the applicant of the present invention has realized a filter circuit 
where the correlation operation at the cell search stage is performed by a high speed operation matched filter, and the 
correlation operation at the path search stage by sliding correlators with a low electric power consumption to cause high 
speed initial synchronization, and reduction in the circuit size and the electric power consumption (Japanese published 

30 unexamined application 19950007-215389). However, a further reduction in the electric power consumption has been 
strongly required. In conventional mobile stations, the path search operation has been performed even at a stationary 
state and a low-speed moving, similarly as at high-speed moving, to cause high electric power consumption. Thereby, 
there have been problems of increase in the device weight caused by use of batteries with large capacity, and insuffi- 
ciency in the battery life. 

35 [0005] The present invention has been made under the above circumstances to offer a DS-CDMA receiver with fur- 
ther reduction in the circuit size and the electric power consumption. 

Disclosure of Invention 

40 [0006] A DS-CDMA receiver according to the present invention has smaller number of multiplication circuits than 
that of the code length of a second spreading code, performs correlation operation at a phase with a center of the cor- 
relation peak train detected at the above matched filter as a center phase, and causes change in the calculation fre- 
quency according to the change in the delay profile. 

[0007] Moreover, a DS-CDMA receiver according to the present invention calculates a correlation by corresponding 
45 number of sliding correlators to that of the range of the phase window with a center of the correlation peak train as a 
center phase to cause change in the width of the phase window according to the width of the correlation peak train dur- 
ing reception. And in the DS-CDMA receiver according to the invention, there is a reduced start frequency of the overall 
path search stage when the change in the delay profile is small at a stationary state and low speed moving, and a higher 
frequency of the overall path search stage when the change is large. 

50 

Brief Description of Drawings 
[0008] 

55 FIG. 1 shows a block diagram of a DS-CDMA receiver of one of embodiments according to the present invention. 

FIG. 2 shows a block diagram of a cell search stage of the embodiment shown in FIG. 1. 
FIG. 3 shows a block diagram of a path search stage of the embodiment shown in FIG. 1 . 

FIG. 4 shows a conceptual diagram of operation of coherent adding at the path search stage of the embodiment 
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shown in FIG. 1. 

FIG. 5 shows a conceptual diagram of search timing at the path search stage of the embodiment shown in FIG. 1. 
FIG. 6 shows a flowchart of a controlling method of a search frequency (a start frequency of the overall path search 
stage) at the path search stage of the embodiment shown in FIG. 1. 
5 FIG. 7 shows a flowchart of a method to decide a width of a phase window of the embodiment shown in FIG. 1 . 

FIG. 8 shows a block diagram of sliding correlators at the cell search stage and the path search stage of the embod- 
iment shown in FIG. 1 . 

FIG. 9 shows a block diagram of a control channel receiving stage of the embodiment shown in FIG. 1. 
FIG. 10 shows a block diagram of a traffic channel receiving stage of the embodiment shown in FIG. 1. 
w FIG. 1 1 shows a block diagram of a transmitting roll-off filter of the embodiment shown in FIG. 1 . 

FIG. 1 2 shows a block diagram of a transmitting stage of the embodiment shown in FIG. 1 . 
FIG. 13 shows a block diagram of a radio frequency control stage of the embodiment shown in FIG. 1. 
FIG. 14 shows a block diagram of a control stage of the embodiment shown in FIG. 1. 
FIG. 15 shows a conceptual diagram of a frame format of a conventional traffic channel. 
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Best Mode for Carrying Out the Invention 



[0009] A DS-CDMA receiver of one of embodiments according to the present invention will be described below, 
referring to the attached drawings. FIG. 1 shows a block diagram of the overall configuration of the above embodiment. 

20 A radio frequency received signal is down-converted to a baseband signal RS. which is further converted to a digital 
signal by an analog digital (A/D) converter 1. The received digital signal is input in parallel to a cell search stage 2, a 
path search stage 3. a control channel receiving stage 4. and a traffic channel receiving stage 6. 
[0010] The cell search stage 2 comprises a matched filter 21 (FIG. 2) to calculate a correlation between the 
received signal and a predetermined spreading code, and detects a synchronous timing with the received signal by the 

25 filter 21 (initial synchronization). The duration to acquire the received signal is called as a chip time, and the matched 
filter 21 completes the correlation operation every one chip time. After detecting the synchronous timing, the stage 2 
identifies a spreading code or a scrambling code peculiar to each base station to select a base station for receiving the 
signals. 

[0011] After the initial synchronization, a path search stage 3 extracts a plurality of predetermined correlation peaks 

30 (for multipass), based on a synchronous timing (frame synchronization and slot synchronization) with the received sig- 
nal, to identify their timings (phases), when multipath signals are received by the multipath phenomenon (a plurality of 
signals with various time differences reach a mobile station by reflection and so on. though one signal is transmitted.). 
A correlator 41 (FIG. 3) comprising a plurality of sliding correlators is used for the above extraction of the correlation 
peaks. The above matched filter has corresponding number of multipliers (not shown) to all the timings required for cor- 

35 relation operation to calculate the correlation for each timing. However, the circuit size and electric power consumption 
required for the above case are large. On the other hand, the circuit size and the electric power consumption of the case 
with sliding correlators are remarkably reduced, comparing with those of the case where corresponding numbers of 
multipliers to all the timings are provided, as much smaller numbers of sliding correlators in comparison with those of 
the latter case, for example, sliding correlators SC1-SCn (FIG. 5) have multipliers ml - mn (FIG. 5). respectively. But. 

40 much time (one symbol duration = spreading factor x chip time) is required for one correlation operation. Then, the cell 
search stage 2 is used only for the initial synchronization, and. after completion of the initial synchronization, the stage 
2 is stopped to reduce the electric power consumption. Furthermore, the cell search stage is started not only for the 
initial synchronization (initial cell search), but also for peripheral cell search. 
■ [0012] The cell search stage 2 selects a correlation peak train with the highest power (for the main path) as the ini- 

45 tial synchronous timing in a subsequent control stage 30 (FIG. 2) after the matched filter 21 . And. the timing is input to 
a path search stage 3 as a control signal CTRLS2. 

[0013] The control channel receiving stage 4 receives control signals for various kinds of controls in the mobile sta- 
tion, despreads the signals by the sliding correlators, using a predetermined spreading code. and. after synchronous 
detection and RAKE combination, extracts control signals from the received signal to demodulate them. And the traffic 
so channel receiving stage 6 despreads the signals by the sliding correlators, using a predetermined spreading code. and. 
after synchronous detection and RAKE combination, extracts traffic channel signals including information signals and 
so on to demodulate them. Multipath signal receiving timings for the control channels and the traffic channels are deter- 
mined based on the search results of the path search stage 3. respectively. 

[0014] The transmitting part of the mobile station is provided with a transmitting stage 8 to generate transmitting 
55 signals obtained through multiplication of the data to be transmitted by the spreading code, and transmit them from a 
radio frequency stage (not shown) after waveform shaping of them by a transmitting roll off filter 7. The radio frequency 
stage is controlled by a control stage RFC to output control signals CS. DT to the radio frequency stage. 
[001 5] After despreading. and RAKE combination of a signal at the control channel receiving stage 4 and the traffic 
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channel receiving stage 6. the signal is input to the receiving buffer 5. and processed by a digital signal processor (DSP) 
through a DSP bus B1 and a DSP interface 1 3. The signal to be transmitted from the above transacting stage 8 is input 
to a transmitting buffer 9 through the bus B1 for temporary storage, and then to the transmitting stage 8 from the buffer 
9. 

5 [0016] The outputs of the cell search staae 2 and the path search stage 3 are processed by a microprocessor (not 
shown) through a microprocessor bus B2. and a microprocessor interface 1 2. A control stage 1 1 is provided for other 
controls in the mobile station, and is similarly driven by the microprocessor (MPU). too. The receiving buffer 5 is also 
connected to the bus B2. and processed by the microprocessor. 

[0017] FIG 2 shows details of a cell search stage. An input signal SI2 (corresponding to the output of the A/D con- 
w verier 1 in FIG 1 ) is input to a matched filter 21 . and a correlator 27 comprising a plurality of sliding correlators. A short 
code register 25 is connected to the matched filter 21 and the correlator 27. Moreover, a long code (scrambling code) 
generator 26 which may be a long code register is also connected to the correlator 27. 

[0018] Identification of each cell, that is. the identification of a base station for receiving signals is performed with 
the long code (scrambling code), and a possible long code is generated by the long code generator 26. However, direct 
15 comparison with all the long codes requires much time for the identification. Then, the long codes are d.v.ded into a plu- 
rality of long code groups to identify the long code groups in the first place. The long code groups are identified w.th 
identifying codes for long code group which are supplied from the short code register 25. However, it is necessary to 
detect a slot synchronization timing before identification of the long code groups, and perform a correlation operation 
between the received signals and a synchronous short code, which is common to all the cells, and is also supplied from 

20 the register 25. eio,„,i 
[0019] The cell search stage 2 sequentially performs a correlation operation between the received signals SI2 and 
a synchronous short code common to all the cells by the matched filter 21 . in the first place, to calculate the powers of 
correlation signals of the outputs of the matched filter 21. and, then, recursive integration of the correlation powers .s 
performed for one slot duration in a recursive integration circuit 22. The recursive integration of the correlation powers 

25 is performed for one or more frames to treat a timing with the highest output correlation power as a slot synchronous 
timing and also as a receiving timing of a main path. According to the timing of the main path, a correlation operation 
between the short codes for identification of the long code groups and the received signals is performed, using a plu- 
rality of sliding correlators, to calculate the powers of the correlation outputs, using a correlation power operator 28. 
Long code groups are identified, based on the correlation powers using a plurality of short codes for identification of the 

30 long code groups. After identifying the long code groups, a long code of a base station (cell) for receiving signals is .den- 
tified. based on the correlation output powers between the received signals and the combined codes of all the long 
codes in the groups with the short code for a symbol duration assigned to the control channel. 

[0020] FIG 3 shows details of the path search stage 3. A received signal of pilot blocks (from PILOTk-1 to 
PILOTk+1 in FIG 4) in an input signal SI3 (corresponding to the output of the A/D converter 1 in FIG. 1) is despread in 

35 a correlator 41 comprising a plurality of sliding correlators after storage in an input buffer 40. Coherent adding of the 
output of the correlator 41 is performed in a coherent adder 42 after elimination of the modulation effect of the pilot sig- 
nals Thereafter the correlation powers of the coherent adding outputs are calculated in a correlation power operator 
45 The correlation powers are averaged with a recursive integration circuit 46. and input to an output buffer 47. A con- 
trol circuit 48 is connected to the output side of the circuit 46. and detects a multipath timing of the output of the circuit 

40 46 The above control signal CTRLS2 is also input to the circuit 48. The correlator 41 comprises smaller number of slid- 
ing correlators than that of the code length of spreading codes for the correlation operation. Timings, during which the 
sliding correlators are applied, depend on timings of the control signal CTRLS2 and the detected multipaths. The oper- 
ation frequency of the correlator 41. that is. the frequency to calculate the correlation is adjusted according to change 
in the above delay profile, wherein the coherent adder 42 comprises averaging circuits 43. 44 in parallel, which average 

-<5 the pilot symbols of the adjacent symbol durations, respectively, to calculate the powers for those averages in the cor- 
relation power operator 45. In addition, the results of the power calculations are averaged at the recursive integration 
circuit 46 to detect the latest multipath timing from the average results. The detected multipath timing is applied to the 
receiving timing of the correlators for receiving the control channel and the traffic channel from the top of the slot just 
after the detection. 

50 [0021] FIG 4 shows a conceptual diagram of an operation of coherent adding at the path search stage 3. An aver- 
age value Avk is generated by averaging pilot symbols of pilot blocks PILOTk-1. and PILOTk in adjacent slots SLOTk- 
1 and SLOTk Similarly, an average value Avk+2 is generated by averaging pilot symbols of pilot blocks PILOTk+1 and 
PILOTk+2. respectively in the subsequent slots SLOTk+1 and SLOTk+2. The correlation powers are calculated from 
the average values to obtain the average value AVPk+2. 

55 [0022] Thereby it may be possible to eliminate noises and interferences and realize improvements in detection 
accuracy and stabilization of multipaths by coherent adding to calculate the averages for a plurality of slots. And the 
operation may be performed with comparatively small number of sliding correlators, as the objects of the coherent add- 
ing are only the pilot symbols, and only despreading of the received signals of the pilot symbols are required. Moreover. 
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as described below, the range of the phases (phase window), and the calculation frequency for calculation of delay pro- 
files may be adjusted to reduce the number and the frequency of usage of sliding correlators to the minimum. Though 
the pilot blocks of a traffic channel shown in FIG. 4 are object received signals for calculation of the delay profiles of the 
path search stage, pilot blocks of control channels, also, may be the object received signal for calculation of the delay 

5 profiles, if the pilot blocks of the control channel have similar configuration. 

[0023] FIG. 6 shows one of flowcharts for a controlling system to adjust the start frequency of the above path 
search stage 3 as a whole. At the step S1. the number of multipath correlation peaks in a predetermined duration (in 
one symbol duration) with a level equal to or exceeding a predetermined level is counted, and it is judged whether the 
multipath number is changed from that of the previous multipath detection result. When the multipath number is 

w changed, an operation frequency F is set to the maximum value Fmax at the step S2 to end the operation. On the other 
hand, when the number is not changed, overall evaluation of the phase shifts of each multipath is performed at the next 
step S3. One example of the evaluations will be described below. It is assumed that the number of multipaths is r, the 
multipaths are put as P1-Pr. the peak levels of each multipath signal as Ps1 - Psr. and the phase shifts as 51 - 5r. Then, 
the" total sum of phase shifts PD. and a value PDR by dividing the PD by the total sum of the peak levels are used as 

?5 indices. The indices are expressed by formulae (1). (2): 

;=1 . 

20 



30 



PDR = 



(-1 



(2) 



When the indices such as PD and PDR are represented by X it is judged at the step S4 whether A is exceeding a pre- 
determined threshold Ath. If not exceeding, the processing is stopped as it is. and if exceeding, the frequency is cor- 
reeled by a frequency correction function AF(A-Ath) to {F+F(A-Ath)} at the step S5. As the frequency correction function. 
35 for example, there may be used a function, by which a frequency change aF is generated through multiplying (A-Ath) 
by a predetermined proportional constant, and the frequency F is changed in a range between the minimum value Fmin 
and the above Fmax. 

[0024] FIG. 5 shows a conceptual diagram of timings during which the above correlator 41 is applied. When the 
above received signals SI3 are expressed by a time sequence signal Slk. the length of each signal is one symbol dura- 

jo tion, and each sliding correlator SC1-SCn calculates the correlation between the signals Slk and a predetermined 
spreading code at sequentially shifted timings. The operation interval CP1-CPn of each sliding correlator SC1-SCn is 
equal to the symbol duration, and the correlation is calculated within a range CRR with a some margin around a center 
CEN (shown by a dotted and dashed line in the signal Sk) of each symbol duration. Thereby, the synchronization is 
recoverable without overlooking peak positions, even when there is a change in signals. 

45 [0025] FIG. 7 shows a flowchart of a method to decide a width of a phase window. Extracted correlation powers are 
sorted after extracting correlation powers equal to or exceeding a predetermined value (Step S71). wherein predeter- 
mined number of correlation peaks are extracted from a correlation power with the maximum level to generate a corre- 
lation peak train where the correlation peaks are arranged in a time sequence (Step S72). Phases Tb. Te at the 
beginning and the end of the correlation peak train are obtained, and a time distance between them AT=Te-Tb is calcu- 

50 lated (Step S73). Thereby, a range for the correlation operation with the sliding correlators, that is. data for decision of 
the phase window may be obtained. 

[0026] The width W of the phase window is obtained through multiplying AT by a predetermined coefficient w (1.0 
or more), and the window is disposed in a range with an equal width in both directions around a center between Tb and 
Te. that is. (Tb+Te)/2 (Step S74). wherein the value of w depends on conditions of the communication channels, and is 
55 set larger under the poor conditions. 

[0027] The above sliding correlator SC1 comprises a multiplier ml . and an integrator sum1 to integrate the output 
of the multiplier. In the multiplier ml. multiplication of a predetermined spreading code and received signals is per- 
formed. Corresponding number of sliding correlators to the maximum width of the phase window (one symbol duration) 
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are provided When the width of the phase window is narrower than that of the maximum width, unused sliding correla- 
tors are stopped to reduce the electric power consumption. As the sliding correlators SC2-SCn have s.mriar configura- 
tions descriptions of them will be eliminated. 

f00281 A block diagram shown in FIG. 8 describes that the correlator 27 of the cell search stage 2 and the correlator 
il of the path search stage 3 make common use of the sliding correlators. As shown in FIG. 8. the number n o sl.d.ng 
correlators Sd-SCn are used for the correlator 41 . and sliding correlators SCI-SCm of them are also used for the cor- 
relator 27 Multiplexers MUX1-MUxm are connected to the output side of the sliding correlators SC1-SCm. respectively. 
The outputs of the multiplexers are connected to either a correlation power operator 28 of the cell search stage 2 or the 
coherent adder 42 of the path search stage 3. Thereby, the cell search stage 2 and the path search stage 3 make com- 
mon use of the sliding correlators SC1-SCm to reduce the overall circuit size remarkably in comparison with the case 
where sliding correlators are provided independently without their common use. 

r00291 FIG 9 shows a block diagram of details of the control channel receiving stage 4 comprising a sliding corre- 
ctor 61 to which a received signal SI4 (corresponding to the output of the above A/D converter 1) is input; a synchro- 
nous detection circuit 63 to perform synchronous detection of data symbols using the results of channel estimat.o f om 
the outputs of the sliding correlator 61 based on the pilot symbols; and a symbol adder 64 to perform RAKE combination 
of the outputs of the synchronous detection circuit 63. A register 62 to supply a predetermined long code and short code 
(spreading code) is connected to the sliding correlator 61 . A long code is supplied to the sliding correlator 61 at a timing 
designated by the above control signal CTRLS3, that is. a timing detected at the path search stage 3. Moreover the 
control signal CTRLS3 is also input to a symbol adder 64 to perform RAKE combination for integration of all. the multi- 

r003oT nal FIG 10 shows a block diagram of details of the traffic channel receiving stage 6 comprising a sliding corre- 
ctor 71 to which a received signal SIS (corresponding to the output of the above AID converter 1 ) is input; a synchro- 
nous detection circuit 74 to perform synchronous detection of data symbols using the results of channel estimation from 
the outputs of the sliding correlator 71 based on the pilot symbols; and a symbol adder 75 to perform RAKE combination 
of the outputs of the synchronous detection circuit 74. Registers 72, 73 respectively to supply a predetermined short 
code (spreading code) and long code (scrambling code) are connected to the sliding correlator 71 Tha spreading 
codes are supplied to the sliding correlator 71 at a timing designated by the above control signal CTRLS3, hat s. a tim- 
ing detected by the path search stage 3. An EX-OR operation of the long code and the short code supplied from the 
registers 72 73 is performed every bit. and the rasults are supplied to the sliding correlator 71 as a combined code^ 
Moreover, the control signal CTRLS3 is also input to a symbol adder 75 to perform RAKE combination for integration of 

m03lV° rr FIG° n i2 shows a block diagram of a transmitting stage 8. The transmitting stage 8 has a spreading code 
modulation circuit 91 to which an input signal SI7 (corresponding to the output of the above transmitting buffer 9) is 
input and registers 92. 93 respectively to supply short codes and long codes are connected to the spreading code mod- 
ulation circuit 91 . Similarly as the receiving side, combined codes of the short and long codes are generated to spread 
the input signals by the combined codes. 

100321 FIG 11 shows a block diagram of a transmitting roll-off filter 7. The transmitt.ng roll off filter 7 having an FIR 
filter 81 eliminates high frequency components of the input signal SI6 (corresponding to the output signal of the trans- 
mitting stage 8). The output of the FIR filter 81 is converted to an analog signal by a D/A converter 82 as an output s,g- 

rOuS] 6 FIG 13 shows a block diagram of a control stage 10 to control a radio frequency control stage. The stage 
10 comprises a control circuit 101 to control a gain of an amplifier for received signals; a control circuit 102 to control a 
gain of an amplifier for transmitting signals, a control circuit 103 for frequency control of a radio signal receiving stage; 
and a converter 104 for general-purpose A/D and D/A conversion. The control circuits 101-103 outputs control signals 
CS1 . CS2. and CS3. respectively. A data signal DT8 is input to the converter 104. and output from the above converter 

ro°034] FIG 14 shows a block diagram of a control stage 11 . The stage 1 1 has an interrupt controller 1 1 1 connected 
to an MPU bus B2; and a control register 1 14 to store various kinds of microcodes. A timer 1 12 and a counter 1 13 are 
connected to the controller 111 and the control register 114 to set control timings. 

Industrial Applicability 

[0035] As mentioned above, the DS-CDMA receiver according to the present invention has an excellent advantage 
to reduce the circuit size and the electric power consumption. 



Claims 



1 . A DS-CDMA receiver comprising: 
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a cell learch stage performing detection of a timing for an initial synchronization with a received signal; and 
a path search stage comprising smaller number of s.iding correlators than the code length of spreading code, 
and capiat ng a delay profile of said received signal after beginning of calculate of a correlat.on w,th a 
Reading code at a pise having an initial synchronization position detected at said cell search stage as a 

characterized in that said path search stage causes change in the calculation frequency or the star, frequency 
of the overall path search stage according to change in said delay profile. 

9 A ns CDMA receiver according to claim 1. characterized in that, in said path search stage, the start frequency of 
2 ' L D oveX s elrch sTage is'educed low in the case of the small change in said delay profile, and the frequency 
of the overall path search stage is controlled high in the case of the large change in sa,d profile. 

3 A DS-CDMA receiver according to claim 1 , characterized in that, in said path search stage, the change in said delay 
profile depends on the number of the correlation peaks of multipath s.gnals. 

peaks of multipath signals by the total sum of the correlation peaks. 
6. A DS-CDMA receiver comprising: 

. r«u ^parrh stace detectinq a timing of an initial synchronization with a received signal; and 

30 said detected multipath widths. 

7 - A DS CDMA receiver according to claim 1. or 6. characterized in that, in said path search stage, an object for cal- 
cuSon oTsJd'rrelation is P L blocks of a control channel or those of a traffic channel, or those of a common 
pilot channel of received signals. 

35 
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FIG. 1 
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FIG. 2 
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FIG. 4 
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FIG. 5 
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FIG. 6 
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FIG. 7 
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Sorting of correlation peaks equal to or 
exceeding a predetermined value in order of size 



Generation of correlation peak train in time 
sequence, after extraction of predetermined 
number of correlation peaks from correlation 
power with maximum level 



Calculation of time distance AT=Te-Tb 
using phases Tb, Te at beginning and end 
of correlation peak train 
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FIG. 10 
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FIG. 1 1 
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FIG. 1 3 
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